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We investigate theoretically the Raman coupling between two internal states of a trapped low- 
density quantum-degenerate Fermi gas. In general, the trap frequencies associated with the two 
internal states can be different, leading to the onset of collapses and revivals in the population 
difference AA'^ of the two internal states. This behavior can be changed drastically by two-body 
collisions. In particular, we show that under appropriate conditions they can suppress the dephasing 
leading to the collapse of AA'^, and restore almost full Rabi oscillations between the two internal 
states. These results are compared and contrasted to those for a quantum-degenerate bosonic gas. 

PACS numbers: 03.75.Fi, 75.45.-|-j, 75.60.Ej 



I. INTRODUCTION 

Quantum-degenerate samples of low-density fermionic 
atomic gases [y, |^ are arguably even more interesting 
than their bosonic counterparts due to the funda- 
mental role played by the Pauli exclusion principle in 
their dynamics. This principle renders their experimen- 
tal realization particularly difficult, since in the simplest 
case, the collisions that are essential in evaporative cool- 
ing largely disappear as the temperature of the sample 
goes to zero. For this reason, more elaborate techniques, 
involving e.g. the use of several isotopes, or sympathetic 
cooling via a bosonic system, have been used to achieve 
degeneracy J . Collisions also leave behind holes in the 
Fermi sea |4 ; these holes are difficult to fill and are be- 
lieved to limit the temperatures that can be achieved to 
about 0.2 Tp, where Tp is the Fermi temperature. 

From a theoretical viewpoint, Fermi systems also 
present a number of difficult challenges. In particu- 
lar, they are not amenable to a mean-field description. 
Hence they cannot be analyzed in a classical-like for- 
malism such as the Gross-Pitaevskii equation, which has 
proven remarkably powerful in describing many aspects 
of quantum-degenerate bosonic systems. On the other 
hand, the additional complexity of Fermi systems also 
offers much promise. One can hope to be able to manip- 
ulate them into strongly non-classical states, with poten- 
tial applications in atom interferometry. Also, the non- 
linear mixing of fermionic matter waves is expected to be 
very different from the bosonic case. As such, Fermi sys- 
tems promise the extension of nonlinear atom optics |^ 
to a regime without counterpart in traditional nonlinear 
optics. 

Just as is the case for bosons, a cornerstone of the 
manipulation of fermionic matter waves is their interac- 
tion with light. In this paper, we discuss the specific 
situation where transitions between two internal states 
of a trapped quantum-degenerate Fermi system at zero 
temperature are induced by Raman coupling. A straight- 
forward generalization of this model could be used to de- 



scribe an output coupler for an atom laser Q|. Our goal 
is two-fold: first, to understand the difference between 
the bosonic and fermionic dynamics; and second, to de- 
termine the role of two-body collisions on the evolution 
of the system. 

We first consider the dynamics of the system in the 
absence of collisions. Section II introduces our model 
and derives the Heisenberg equations of motion for the 
relevant atomic fields, and section III compares the re- 
sulting dynamics with those for a corresponding Bose 
gas. Collisions are introduced in section IV. The result- 
ing equations of motion are solved numerically in the 
framework of a time-dependent Hartree-Fock theory for 
the case of fermions, and in a standard mean-field theory 
for a bosonic sample. Again, we give a detailed compari- 
son of the two situations, and illustrate how collisions can 
change the fermionic dynamics in a non-trivial fashion. 
Finally, Section ^ is a summary and conclusion. 

II. MODEL 

We consider a two-component quantum-degenerate 
atomic system trapped in a one-dimensional, harmonic 
potential with each component corresponding, e.g., to 
one internal hyperfine spin state. In general, the cou- 
pling of the atoms to the trapping field is different for 
the two (spin) components |-|-) and |— ), so that they see 
trapping potentials of different frequencies lu^ and w_. 
The two internal states are coupled by a Raman-type in- 
teraction of frequency v equal to the spin-flip transition 
frequency of the atoms in the ground state of the two 
trapping potentials. This model, which is summarized 
on the diagram of Fig. |l], is described by the second- 
quantized Hamiltonian 



+ hg [ dx \e-''''i/l{x)^_{x) + h.c] , (1) 
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FIG. 1; Two-component Fermi gas in a harmonic trap. The 
trapping frequencies and lo- correspond to the two in- 
ternal states, which are coupled via a spin-flip transition at 
frequency v resonant with the frequency difference of the trap 
ground states. 



where % is the Raman coupling strength. The first- 
quantized Haniiltonian describing the trapping potentials 
associated with the internal states 1-1-) and I—) is 
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2m dx^ 

with E± being the energy of the internal state |±) 
Raman resonance condition is therefore, with Ua = 

E-, 



(2) 
. The 



(3) 



We remark that this resonance condition effects only 
the coupling between the two trap ground states: the 
coupling between all other levels is off-resonant for u)+ ^ 
Lo^. Hence, introducing a small detuning even for the 
ground states does not significantly alter the dynamics 
of the system. 

The atomic field operators corresponding to the 
two traps obey the fermionic, respectively bosonic 
(anti)commutation relations 



^{x),^{x') 



5ij5{x - x'), 



= 0, 



(4) 



where i,j = { + , — }. 

For the harmonic potentials at hand, the Heisenberg 
equations of motion for the atomic field operators take 
the same form, independently of whether the atoms are 
bosonic or fermionic. It is convenient to expand them in 
terms of eigcnstates {un{x)} of one of the trap Hamilto- 
nians H±, say, iJ+ for concreteness, as 



^+ix,t) 



E 



Un{x)an{t), 



(5) 



where the a„'s and 6„'s satisfy either fermionic or bosonic 
commutation relations. In both cases, this expansion 
readily yields the Heisenberg equations of motion 



I—— = AnUn + gbn, 

dr 

i-T- = Bnbn + Cnbn+2 + Dnhn-2 + fffln, 

dr 

where we have introduced the coefficients 
1 
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and the ratio 
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\ {P^ - 1) v/(n + 2)(n+l). 



- - I) ^n{n - I), 



(3 = LU+luJ- 



(7) 



(8) 



of the trap frequencies. The dimensionless time r is 
scaled to r = w+t, and so is the dimensionless cou- 
pling strength g — g/ijJ+. 

We emphasize that while the operator a„ describes the 
annihilation of atoms in level n of the upper trap, a sim- 
ilar interpretation of the 6„'s is not possible, since they 
result from the expansion of the field operator of atoms 
in the internal state |— ) on the basis of the "-|-"-trap. 
Denoting the eigcnstates of the single-atom Hamiltonian 
of the lower trap as {vnix)}, the "true" annihilation op- 
erators Cn associated with the trapped atoms in the |— ) 
internal state are related to the &„'s by the mapping 



(9) 



where the mapping matrix element Tnm is the overlap 
integral 



T — 



dxVn{x)u„i{x). 



(10) 



III. DYNAMICS 



In this section, we compare the dynamics of ideal non- 
interacting bosonic and fermionic systems evolving under 
the influence of the Raman coupling. We proceed by nu- 
merically solving the Heisenberg equations of motion (||) 
for a sample of N atoms initially in the internal state | +) 
and at temperature T = 0. For bosonic atoms, all atoms 
are therefore initially in the "-|-"-trap ground state, while 
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for fermions they fill the lowest N trap levels. The cor- 
responding initial states are correspondingly 

1=0 

in the case of fermions, and 

|VB(0)) = ^ar|0)+®|0)_ (12) 

for bosonic atoms. 

We consider, first, the case of noninteracting fermionic 
atoms. For trap frequencies approximately equal, /3 ~ 1, 
Eq. (^ suggests the existence of two limiting situa- 
tions, at least in the case of fermions. (We will revisit 
this point when discussing low-temperature bosonic sys- 
tems.) In the first one, which we call the "strong-coupling 
regime" in the following, g « TV, so that the inter-trap 
coupling dominates the dynamics and the intra-trap cou- 
pling terms bn±2 can largely be ignored. In contrast, the 
"weak-coupling regime" g <^ N is dominated by intra- 
trap coupling. 

As a first measure of the system dynamics. Fig. ^ 
shows the difference 

AN{r) = ^ j dx [(*;(.T)*+(x))-(*L(a;)*_(.T))" 
- ^ E ((4 Wa„(r)) - {ci{r)cn{T))) . (13) 

n 

between the populations of the "-I-" and "-" traps. Fig. 
||a is for the strong-coupling regime, and Fig. ^3 for the 
weak-coupling regime. 

One can gain some intuitive understanding of the 
strong-coupling regime by remarking that in that regime, 
intra-trap transitions remain small, so that the Raman 
coupling is predominantly between levels of the two traps 
with equal quantum number n. To lowest order, these 
transitions are all at the Rabi frequency g. However, this 
simplest description cannot explain the result of Fig. ||a. 
Rather, it is necessary to include at least their lowest- 
order corrections, i.e, 

-\j~9' + \ {An B^f c g + (^^^) n'. (14) 

Such an n-dependence of Rabi frequencies is known to 
lead to collapse and revival phenomena, as was first dis- 
cussed in the context of the Jaynes-Cummings model [Q, 
where fin oc ^/n. This is precisely the type of behavior 
exhibited by AiV in the strong-coupling regime. Because 
of the n^-dependence of ri„, it is expected that the lowest 
trap levels, i.e. the atoms in the deep Fermi sea, play a 
dominant role in the appearance of the revivals. We ver- 
ified that the populations of the lowest trap levels indeed 
oscillate more or less in phase, while those of higher n 
levels dephase rapidly. 
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FIG. 2: AN{t) for N = 10 fermions and trap ratio f3 = 0.9, as 
obtained from a numerical integration of Eqs. (^): (a) strong 
coupling regime with g = 10.0; (b) weak coupling regime with 
g = 1.0. The dimensionless time r is in units of 1/uj+. 



We remark that both collapses and revivals of AA^ dis- 
appear when the two trap frequencies are identical, since 
for P — 1, we have An — Bn = n and hence Qn = 9- In 
addition, intra-trap transitions vanish in that case, due 

to Cn = D^ = 0. 

Fig. 1^ shows the inversion AiV between the total trap 
populations in the weak-coupling regime, g <ti N . In this 
limit, inter-trap and intra-trap coupling occur on similar 
timescales. Immediately following a Raman transition 
from the |+) to the |— ) internal state, the population of 
the "-"-trap starts to undergo a redistribution between 
its levels. The combined effects of the intra- and inter- 
trap transitions result in that case in a random-looking 
evolution of AN{t) of Fig. |b. 

We now briefly turn to the case of a Bose gas. For a 
sample at zero-temperature, T = 0, and initially in the 
internal state |-|-), all atoms are in the ground state of 
the "+"-trap at t = 0. As a result, the strong-coupling 
regime is characterized by almost perfect Rabi oscilla- 
tions of the atomic population between the two trap 
ground states, with a very small fraction of the atoms 
coupling to higher modes due to intra-trap transitions. 
This behavior is also largely preserved for — 1 ^ 
g ^ N. (We recall that for fermions the right-hand 
side of this inequality corresponds to the weak-coupling 
regime, dominated by intra-trap transitions.) This differ- 
ence between bosons and fermions can readily understood 
from Eqs. (^, which show that in the T = bosonic 
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case, intra-trap coupling first occurs between the levels 
n = and n = 2 of the "-"-trap, with coupling coefficient 
Z?2 = [0^ — l)/4- As long as this coupling remains small 
compared to the inter-trap coupling g, the system acts 
effectively as a two-mode system. In other words, for low 
temperature bosonic systems, the weak-coupling regime 
is not characterized by 5 w as is the case for fermions, 
but rather by g ^ — 1. 

As is to be expected, the difference between fermionic 
and bosonic systems is reduced as T is increased. At 
first, the sharp edge of the Fermi-Dirac distribution soft- 
ens, resulting in slightly reduced (strong-coupling) col- 
lapses and revivals of the fermionic system. On the other 
hand, for T 7^ 0, bosons occupy higher trap states, re- 
sulting in a spread in Rabi frequencies participating in 
the population difference signal. This in turn leads to 
collapses and revivals rather than the perfect T = Rabi 
oscillations. Increasing the temperature further leads of 
course to undistinguishable behaviors of the bosonic and 
fermionic systems. 



IV. COLLISIONS 

In this section, we discuss the effect of collisions on the 
preceding results. Collisions are of course central to the 
dynamics of quantum-degenerate atomic systems. They 
are essential in the evaporative cooling of the sample, and 
also provide a nonlinearity that can lead to the nonlin- 
ear mixing of matter waves. In bosonic systems, much 
new physics can be studied, e.g. by changing the sign of 
the scattering length of s-wave collisions. In the case of 
fermions the creation of holes in the Fermi sea, and the 
filling of these holes by additional collisions results in a 
heating that appears to fundamentally limit the temper- 
atures at which these samples can be cooled First, 
we discuss the way collisions impact the operation of the 
Raman coupler in the case of fermions, and later compare 
these results with those for a bosonic sample. 

It is well known that in fermionic atoms, the Pauli 
exclusion principle forbids the existence of s-wave scat- 
tering between atoms in the same internal state. In ad- 
dition, p-wave scattering is generally negligible. Hence 
two-body collisions are described by the Hamiltonian 

-Hcoi^Uo J dx^l{x)i'l{x)^^{x)^+{x), (15) 

where Uq — Anh'^ap/m is the interaction strength with a 
being the s-wave scattering length and p the character- 
istic density of the system. Again, we expand the field 
operators according to Eq. (^ in terms of the basis 
and obtain 



i.j,k,l 

where the matrix element 



ijki a\b^jbkau 



Uijki = / dxui{x)uj{x)uk{x)ui{x) 



(16) 



(17) 



characterizes the scattering between different levels. We 
note that Uij^i is symmetric under permutations. 

In the presence of this quartic Hamiltonian, the Heisen- 
berg equations of motion for the operators a„ and 6„ in- 
volve cubic combinations of operators. To close this sys- 
tem of equations, we invoke a time-dependent Hartree- 
Fock ansatz, which has proved to be very successful 
in the treatment of many-particle quantum systems 
to factorize products of operators, of the generic form 



h]{T)b,(T)ak{T), hy 



bl{T)b,{T)au{T) « (5l(T)6,(r))a,(T) - (6j(T)afe(r))6,(r), 

(18) 

where the expectation value is over the state \iI)f{^)) 
since we work in the Heisenberg picture. At this level 
of approximation, we neglect all contributions from pair- 
ing. This factorization scheme readily yields the time- 
dependent Hartree-Fock equations of motion (in dimen- 
sionless variables) 



.dan 
.dbn 



[{AnSnk + Q"!) dk ~ {Qtk* - ~g5nk) h 

k 

^ [{BnSnk + Qnk) h ^ {Qnl " 9^nk) Ofc 



+ Cnbn+2 + Dnbn-2, 



(19) 



where we have introduced the time-dependent coefficients 

Qnlir) = UoJ2UmM^Kr)a,{T)) 

hi 

Qnkir) = UoJ2^mMbMb,{T)) 

Qnkir) = UoJ2Un,,k{dKT)b,{T)), (20) 



and Uq = Uo/huj^. is a dimensionless interaction strength. 

The effect of collisions is illustrated in Figs. ^,b, 
which show the population inversion AN{t) for two val- 
ues of the interaction strength Uq. 

For weak enough collisions, the dynamics of the sys- 
tem is not significantly altered, as should of course be 
expected. However, we observe a quantitative change in 
the dynamics of AA^ as C/q is increased. Instead of a 
collapse and revivals, AN{t) now undergoes nearly full 
Rabi oscillations. 

A first hint at the cause of this changed behavior is 
offered by Fig. U, which shows a snapshot of the level 
opulations in the "-|-"-trap for the cases of Fig. ||a and 
respectively. We observe that the smaller value of Uq 
corresponds to an inhomogcncous level population distri- 
bution, whereas the higher nonlinearity causes the trap 
levels to be almost equally populated. 

A more quantitative understanding of the role of col- 
lisions can be gained by estimating how the nonlinear 
terms in Eq. ( |l^ ) modify the (coUisionless) Rabi fre- 
quency. A numerical evaluation of the coefficients Uijki 
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FIG. 3: Population difference AA'^(r) for A'^ = 10 fermions, a 
trap ratio /3 = 0.9 and in the strong-coupling regime g — 10.0. 
The plots, which result from the numerical integration of Eqs. 
(^^, are for different strengths of the two-body collisions: (a) 
Uo ~ 0.01; (b) Uq = 0.1. Dimensionless time t in units of 

1/L0+. 



shows that elastic collisions, i ^ j = k = I, dominate 
the dynamics of the system. In addition, U„nnn turns 
out to be a decreasing function of n. Keeping the elas- 
tic contributions to the collision-induced dynamics only, 
and neglecting as in the strong-coupling regime of section 
III the effects of intra-trap coupling terms 6n±2, one find 
that as a result of collisions Eq. (|lj) is approximately 
changed to 



(21) 

Fig. H shows, as a function of Uq, the time-dependent 
Rabi frequencies fi^^ (r) averaged over a time interval Q 
large compared to their inverse, 

^n = ^ j'^ drn^^^T). (22) 

Because is a decreasing function of n, its contri- 

bution tends to compensate the ri^ dependence of Eq. 
(|l^). As a result, there is a range of collision strengths 
for which the dependence of VL^^{t) on n largely dis- 
appears. In this range, paradoxically, the dynamics of 
the collision-dominated Fermi system resembles that of 
a coUisionless Bose system. From this admittedly crude 
argument - which is however consistent with our full nu- 
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FIG. 4: Occupation of the first 15 upper trap levels at the 
dimensionless time r = 490/a;+, (a) for the parameters of 
Fig. ^a; and (b) for the parameters of Fig. ^d. 
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FIG. 5: Semi-logarithmic plot of the time-averaged gen- 
eralized Rabi frequencies for the six lowest trap levels as a 
function of the nonlinear parameter, otherwise the same pa- 
rameters as in Fig. ^ are used. Q,„ is given in units uj+ , JJo in 
units of TujJ+. 



merical results - we also conjecture that for even larger 
C/q, the approximate cancellation of the rt-dependence of 
the Rabi frequencies will disappear and we expect an 
overall dephasing and decay of the population difference 
AiV(i). It has unfortunately proven prohibitive to try 
and check this conjecture numerically. 

We now turn to the case of bosonic atoms. Bose statis- 
tics allows for s-wave collisions between atoms in the 
same spin state, so that the collisional Hamiltonian is 
now 

Hcoi = u+ [ dxiil{x)i>\{x)i>+{x)i>+{x) 
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+ 2t,./..4U.)*L(x)*,w*_w, (23, 

where the Ui, i — {+,— ,a;} characterize the strength of 
the colhsions. In the foUowing we assume for simphcity 
J7+ = U- = Uq and Ux — rj^Uo. 

To truncate the Heisenberg equations of motion for 
the field operators, we now invoke a mean- field approx- 
imation, factorize all products of operators, and replace 
the resulting expectation values by time-dependent c- 
numbers. This gives 



. d{a„) 

' dT 



where 



k 

- J2 (BnSnk + Qi'k + v.QZ) ih) + ~g{dn) 



+ Cn{bn+2) + Dn{bn-2), 



(24) 



QTkir) = 2;7o5^f/„,.,,:fc(a,(r))*(a,(T)) 
QI\{t) = 2UoJ2UmMr)r{b,{T)), (25) 

and the expectation values are with respect to the state 
IV'b(O)). Figure H shows the inversion AN{t) for a sam- 
ple of bosonic atoms initially in the internal state |-|-) 
for a nonlinear parameter Uq = 0.5. In contrast to the 
case where collisions are absent and we have full Rabi 
oscillations, c.f. Sec. Ill, here one starts observing a 
damping of the oscillations. This is clearly a result of 
the scattering of atoms into higher trap states. This is 
illustrated in Fig. |^, which shows the population of the 
first upper trap levels at a fixed time. The transitions 
between the populated trapped states are characterized 
by n-dependent Rabi-frequencies, leading to the onset of 
a dephasing process resembling the situation for nonin- 
teracting fermions We see, then, that in the case of 
interacting bosons, intra-trap scattering is an important 
element of the dynamics of the Raman coupler, which 
rapidly evolves to a multimode behavior; in contrast in 
the intrinsically multimode fermionic case Uijki tends to 
reduce the spread in Rabi frequencies and thus inhibits 
dephasing. 

A remarkable property of the bosonic trap population 
distribution is that only even trap levels are occupied, 
see Fig. |. This is a combined result of three facts: 
(a) at T = all atoms are initially in the ground state 
of one of the traps; (b) s-wave scattering only couples 
trap states of same parity, as expressed by the symmetry 
properties of the collision matrix from Eq. |l^; (c) intra- 
trap coupling only couples trap levels with An ~ 2, see 

Eq. ®. . . _ 

It is known from nonlinear optics |lO[ and atom optics 
Jll| that systems governed by a pair of coupled nonlinear 
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FIG. 6: AiV(r) for = 20 bosons, trap ratio /3 = 0.9, linear 
coup ling g = 7.0 and rjx = 1.5, from the numerical solution of 
Eq. (p4) with Uo = O.S.Dimensionless time r in units of l/uj+ 
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FIG. 7: Level population of the "-|-"-trap at the dimensional 
time r = 25/lj+ for the parameters of Fig. ^ 



Schrodinger equations can reach a regime where the non- 
linear phase shifts dominate their dynamics. Such two- 
mode systems exhibit Rabi oscillations for small nonlin- 
earities, but mode-coupling is inhibited above a certain 
strength of the nonlinearity. This effect is absent in the 
present multimode system, a result of the strong inter- 
mode scattering. 



V. SUMMARY 

The Raman coupling between two internal states of a 
trapped Fermi gas exhibits a rich dynamics, quite differ- 
ent from its bosonic counterpart. This is of course due 
primarily to the fact that a Fermi gas occupies a large 
number of trap states, and hence can never be approx- 
imated as a two-mode system. In particular, we have 
identified two limiting regimes, dependent upon whether 
inter-trap or intra-trap dynamics is dominant. In the 
general situation where the traps associated with the 
two internal states have different frequencies, the first 
of these regimes leads to dynamics characterized by col- 
lapses and revivals. However, two-body collisions can 
under appropriate conditions inhibit this behavior, mak- 
ing the collision-dominated Fermi system more similar to 
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a collisionless Bose system. 

The numerical analysis of fermionic systems appears to 
be presently limited to very small numbers N of atoms, a 
result of the large memory requirements associated with 
the need to keep track of a large number of quantum 
states. Indeed, most of our numerical results are limited 
to N on the order of 10 (especially for the coUisional cal- 
culations), and even this required rather large computing 
facilities. Despite this limitation and its associated lack 
of quantitative predictions, our analysis sheds useful light 
on the dynamics of trapped Fermi systems, in particular 
in the presence of collisions, and will provide useful guid- 
ance in understanding more realistic trapped Fermi gases 
in three dimensions and with a large number of fermionic 



atoms. 
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